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Abstract. An array of superconducting islands with semiconducting nanowires in 
the right regime provides a macroscopic implementation of Kitaev's toy model for 
Majorana wires. We show that a capacitive coupling between adjacent islands leads 
to an effective interaction between the Majorana modes. We demonstrate that 
even though strong repulsive interaction eventually drive the system into a Mott 
insulating state the competition between the (trivial) band-insulator and the (trivial) 
Mott insulator leads to an interjacent topological insulating state for arbitrary strong 
interactions. 
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1. Introduction 

Majorana zero modes are fermions which are their own antiparticles. They are 
believed to exist as effective particles in the middle of the gap of a topological 
superconductor. The recent interest in Majorana zero modes originates in their non- 
Abelian exchange statistics [lj , which is the basis for potential applications in quantum 
computation [2]. Based on the theoretical proposal [3] to realize this exotic states 
in semiconducting nanowires with strong spin-orbit coupling in a magnetic field and 
in proximity to a conventional (nontopological) superconductor, recent experimental 
progress has shown signatures of Majorana zero modes in the tunnelling conductance 
of normal conducting-superconducting [I] and superconducting-normal conducting- 
superconducting systems [5]. More information and references on the fast developing 
field can be found in the recent reviews [6]. 

As the nanowire is one-dimensional interaction effects become important. On the 
one hand, employing field theoretical methods Gangadharaiah et al. [7] have argued 
that strong electron-electron interactions generically destroy the topological phase by 
suppressing the superconducting gap. On the other hand, using a combination of 
analytical and numerical methods Stoudenmire et al. [8] have shown that repulsive 
interactions significantly decrease the required Zeeman energy and increase the 
parameter range for which the topological phase exists. It is believed that the origin 
of these effects is an interaction driven renormalization of the Zeeman gap [U [9]. 
Furthermore, it has been shown that in helical liquids the scattering processes between 
the constituent fermion bands open gaps which in turn lead to a stabilization of the 
Majorana states against interactions [10] and that (an odd number of) Majorana zero 
modes are in fact stable against general interactions [TT]. For further studies of the effect 
of electron-electron interactions on Majorana zero modes in nanowires and two-chain 
ladders see [12] . 

Recently, a macroscopic version of the Kitaev chain (13j, a toy model for a nanowire 
supporting noninteracting Majorana modes, has been proposed [H] in a one-dimensional 
(ID) array of topological superconducting islands. Its advantage over microscopic 
implementations is that the individual parameters of the effective Hamiltonian can 
potentially be tuned in situ. Here, we show that additional capacitances between 
adjacent islands lead to an effective interaction between the low-energy Majorana degrees 
of freedom. We present the phase diagram of the system which demonstrates that 
sufficiently strong repulsive interactions will drive the system from the topologically 
trivial phase into the topological phase supporting Majorana zero modes before 
eventually leading to a Mott insulating state. We discuss how the parameters of the 
system can be tuned by changing the gate voltages and exploit this to propose the 
detection of the different phases and phase boundaries in a tunnelling experiment. 
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Figure 1. A ID array of superconducting islands (light grey) coupled via strong 
Josephson junctions (Ej) to a common ground superconductor. Each island contains 
a pair of Majorana zero modes (white dots) at the end points of a semiconducting 
nanowire (dark grey). The tunnel coupling of individual electrons between the 
superconducting islands is proportional to the energy scale Em- A common gate 
voltage V g can be used to tune the relative strength of the different terms in the 
Hamiltonian. The capacitive couplings between the elements are denoted by C, Cj, 
and C g , respectively. 



2. Model 



We discuss a system consisting of a ID array of N superconducting islands (see figured]). 
Because of the proximity-coupled semiconducting nanowire each of the islands has two 
midgap Andreev states, i.e., Majorana modes, located at the ends of the nanowire [3]. 
We will denote with 7fc a and 7^ the two Majorana operators on island k associated with 
these zero modes. The Majorana operators are Hermitian 7^ = 7^ and fulfil the Clifford 
algebra {7^, 7^} = 25^ u . The total Lagrangian L = T — Vj — Vm of the system consists of 
three terms. Coupling of the Majorana modes on nearby islands leads to the term [T3lfT5] 
V M = i E M J2 k =ilkblk+i a cos[((f) k+1 - <f> k )/2] where <p k is the superconducting phase of 
the k-th island. Here and in the following, we assume for convenience periodic boundary 
condition such that islands 1 and iV+1 are equivalent. Apart from the Majorana modes, 
the term discussed above has the additional degrees of freedom (fi k due to the condensate 
of Cooper pairs. Similar to [H] , we eliminate these by connecting each superconducting 
island with a strong Josephson junction to a common (ground) superconductor. This 
fixes the superconducting phases (up to some quantum phase-slips discussed below) and 
is described by the Hamiltonian Vj = Ej^2 h=1 (l — cos0fc); here, Ej = hl c /2e is the 
effective Josephson coupling of each of the Josephson junctions with critical current J c . 
In the limit Ej ^> Em the junctions effectively pin the phases of all superconducting 
islands to a common value 4>k = 0. As a result Vm reduces to the pure Majorana 
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coupling 

N 

E M = iE M ^2 7fc&7fc+ia, (Ej > E M ). (1) 

k=l 

Finally, the kinetic term T = (h 2 /8e 2 ) £f =1 [C(4 + i ~ 4) 2 + C G fa] + (7i/2e) £f=i # fe 
occurs due to the capacitive couplings, where C denotes the capacitance between 
neighbouring islands and Cq = C g + Cj the (total) capacitance to the ground. Here, 
Cj is the capacitance of the strong Josephson junction and C g the capacitance to a 
common back gate at voltage V g with respect to the ground superconductor. Apart 
from the charging energy, the back-gate introduces a term proportional to the induced 
charge q = C g V g which is tunable with single-electron precision [16] via V g and whose 
effect will become important later. The typical capacitive energy scale Ec = e 2 /2Cs of 
a single islands depends on the total capacitance Cs = 2C + Cg- 

3. Mapping on an effective spin model 

As we have seen the strong coupling to the ground superconductor pins the 
superconducting phase differences to 4>k = and thus changes the energy due to the 
Majorana modes from Vm to Em- The effect of the strong coupling to the ground 
superconductor on the charging energy T is more subtle. We first present the results 
for C = before extending them to nonzero C: as charge and phase are conjugate 
variables the pinning of the phases <pk greatly reduces the effect of charging [13]. An 
effective charging energy still arises due to quantum phase slips through the junctions. 
For example, changing from to 2tt leads to a charging energy (Ej ^> Ec) [17] 

H C k = r A cos[7r(g/e + n fe )] = T A cos(irq/e)V k . (2) 

Here, the tunnelling amplitude is given by Ta — E^ Eje~ SA with Sa = TT 1 Jdr Le = 
a/ 8E j I Ec the dimensionless Euclidean action along the classical trajectory 7a: 0^ G 
[0, 27r] in the inverted potential with Le{t) = —L{t = The cosine term in (T2]) 

occurs due to the Aharonov-Casher interference of the two tunnelling paths <ftk = — > 2-7T 
and <frk = — > — 2tt which lead to an indistinguishable final state and thus interfere 
with a phase difference depending on the total induced charge q + enk [EJ IE] • The 
term = \{\ — Vk) £ {0, 1} is the contribution to the charge due to the parity 
Vh = i'-fkalkb £ { — 1? 1} of the number of electrons on the superconducting island 
encoded in the state of the Majorana zero modes [H2, [T3]. Equation (j2J) is a chemical 
potential term: For large Ta, all the fermionic states are either filled or empty and the 
system is in a band insulating state. 

Going away from this special point and introducing a finite cross-capacitance 
parametrized by rj = 2C/Cy, with rj e [0,1], the classical path for a phase slip in 
<pk does not only involve <pk but also the other phases. To lowest nonvanishing order 



\ Because Vj 3> Vm , we take only the potential Vj into account when calculating the tunnelling action 
S E - 
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in 77, we only need to take into account the change in 4>k-i and 4>k+i and obtain = 
y^8Ej/Ec[l + (7r 2 — 12)r/ 2 /96], which is accurate to more than 2 digits all the way up 
to T] = 1 as numerics confirms. Additionally, the effective capacitance coupling between 
two islands becomes important. This term is generated by simultaneous phase slips of 
the two phases (pk and 4>k+i of neighbouring islands. Due to the symmetry between 
island k and k + 1, we have (pk = (fik+i for the classical path ju : <fik = 4>k+i £ [0,27r] 
with the corresponding Euclidean action Sjj = a/16(2 — rj)Ej/Ec- The Aharonov- 
Casher interference in this case leads to the tunnelling amplitude ~ E^Ej^e -317 
depending on the total charge 2q + e(n& + rifc+i) of the two islands involved. Thus we 
obtain the interaction term 

Huk = IV cos[27rg/ e + 7r(n fc + n fc+ i)] = IV cos(27rg/ e)V k V k +i- (3) 

We stress that OS]) is an interaction term involving four Majorana operators and that the 
amplitude is modulated by twice the induced charge q compared to ©. We call U > 
repulsive interaction as it prefers having an occupied site next to an empty one. For 
strong repulsive interactions the system is driven into a Mott insulator or equivalently 
commensurate charge density wave (CDW) state. 

The complete low-energy Hamiltonian in the limit Ej 3> Em, Eq is given by 
-f^ANNNi = ^2k=i(Hck + H Uk ) + H M , which constitutes the transverse axial next-nearest- 
neighbour Ising (ANNNI) model [19J as can be seen by performing a Jordan- Wigner 
transformation Vk = ilkalkb = erf, ilkblk+ia = a k a t-\-\i resu hing in the spin Hamiltonian 

N 

H ANNNl = ^(Ao-fc + Ua z k a z k+1 + E M ala x k+l ). (4) 

k=l 

Here, the o~ k ,y ' z are Pauli matrices, A = rACOs(7rg/e) and U = IV cos(27rg/e). We 
note that the energies A and U can be tuned via the charge q induced on the 
superconducting islands through the gate voltage V g . As the fermionic Hamiltonian 
without the interaction term proportional to U has been intensively studied before [13], 
we focus here on the effects of the interaction term. In our set-up, this term is most 
important in the case C Cq where 77 fa 1 and Tjj fa exp [— 1 . 25 v/ E j / Eq\ , which 
can be as large as 0.3 for Ej fa Ec; note that in this regime the actions Sa and Sjj are 
still much larger than one such that the semiclassical approximation employed above is 
valid. 

4. Phase diagram 

We first note that the system (TjO) is invariant under Em — > —Em as well as A — > —A 
due to the transformations a k ,y — >■ (— l) k o~ k ' y and a k ,z — > —<J x k z respectively. Thus 
without loss of generality, we assume in the following Em, A > 0. The phase diagram 
of the spin model (T4J) contains four phasesj§| A paramagnetic phase (PM) with a unique 

§ The Hamiltonian Q is brought to standard form by performing the duality transformation fif — 
^lj<k a j> = <J k a k+i- The phase diagram of the resulting model including expressions for the phase 
transitions has been worked out in [19l [22] in the parameters k = U/A and Em/ 'A. 
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Figure 2. Phase diagram of the Josephson junction array The phase diagram 
is invariant under A — > —A. The topological phase is characterized by a doubly 
degenerate ground state. In figure |3] we show the tunnelling conductance along the 
paths ci and c%- Inset: Gap At between the (nearly) degenerate ground states 
(S < 0.02 _Z?m) and the first excited state along the blue dashed line as determined 
by exact diagonalization. 



ground state with (a^) < 0, (af.) = 0. An antiferromagnetic phase (AFM) with doubly 
degenerate ground state and = 0, (of) oc (— l) fe . An "anti phase" (AP) with a 
doubly degenerate ground state with (o%) oc (— l) k and a "floating phase" (FP) between 
the AFM and the AP. For A = 0, the duality transform of the model (J3J) is a sum of two 
quantum Ising chains, while for the noninteracting case it reduces to a single quantum 
Ising chain. 

As we have shown above in our realization of the ANNNI model the parameters 
A and U can easily be tuned via a gate voltage. On the other hand, the coupling 
Em is determined by the overlaps of the Majorana wave functions on neighbouring 
islands and thus fixed by the geometry of the array. Hence it is natural to consider 
the phase diagram as function of A/E M and U/E M , which is shown in figure |2j 
For fixed A/ Em > 1 and sufficiently small interaction U the system is in a trivial 
(band-insulating) phase corresponding to the PM in the effective spin model, which is 
characterized by an unique ground state with (Vk) < 0. By increasing U we cross into 
a topological phase (corresponding to the AFM) with (Vk) = and two degenerate 
ground states |±), distinguished by the (total) fermion parity llfc^-l 1 ' 1 ) = ^l^)- The 
parity protection (also called topological protection) originates from the fact that any 
fermionic perturbation conserves the fermion parity and thus cannot mix the states |±). 
The phase transition between the trivial and topological phase is, for U > 0, located at 

-, 2U = E M UE 2 M 

A A 2A 2 (A-f/)' {) 

In particular we find A = E M + 3U/2 for \U\ < E M and A = 2U for U > E M . 

At large positive U we eventually enter incommensurate and commensurate CDW 
(Mott insulator) states corresponding to the FP and the AP, respectivelylJJ] This region 



| We note that the existence of the FP, and thus the incommensurate CDW state, at small A/ Em has 
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of the phase diagram cannot be reached as long as the induced charge q on the islands 
is homogeneous. However, replacing the common back gate by individual gates for each 
island yields the model (pEJ) with site-dependent parameters Aj = cos(7rgj/e) and 
Ui = IV cos(7r(gj + g i+1 )/e) where = G 9 V^ denotes the induced charge on island i. 
Now using qi = (— l) l q one can enter the Mott phase for q — > e/2. Note that the two 
degenerate ground states in the Mott insulator phase have the same fermion parity and 
thus are not parity protected. 

The main feature of the phase diagram is its strong anisotropy under U — >■ —U. 
While negative interactions suppress the topological phase, for U > the ordering 
tendencies of the first and second term in (jlj) compete with each other. In particular, 
starting from a noninteracting point in the trivial phase, i.e., U — and A > Em, 
competition between the band- and the Mott-insulator will drive the system into the 
topological phase irrespective of the value of A/Em- 

As discussed above the topological phase is characterized by the existence of a 
doubly degenerate ground state with different fermion parity. In a finite system this 
degeneracy is lifted and the value of the resulting gap 5 depends on the number of 
islands N as well as the system parameters A/ Em and U/Em- On the other hand, 
the existence of Majorana end modes is protected by the gap between the (nearly) 
degenerate ground states and the first excited state, which we denote by A T . This 
gap is given by A^ = 4Em at A = U = and decreases when going away from this 
point. However, as we show in the inset in figure [2, it remains of the order of Em 
and significantly larger than 5 for U ^ as long as one stays away from the phase 
boundaries. 

5. Tunnelling conductance 

After presenting the phase diagram, we now turn to its experimental signatures in the 
tunnelling conductance. The parameters A and U can be directly tuned through the 
induced charge q on the islands via a common back gate voltage V g making it possible 
to choose the parameters such that the path crosses one or more phase boundaries. 

Specifically, we consider the path Ci(q) = (IV cos(27rg/e), T/\ cos(iiq/e)) with 
IV = 0.3 Ta = 0.54 Em such that the starting and end points at q = and q = e/2 
lie in the topological phase while the path enters the trivial phase in between (see 
the red curve in figure [2]). In the following we consider an open chain of N islands. 
Following [20], we couple the system to an electronic lead with the tunnel Hamiltonian 
H T = t'-fi a J2ka( c k(T ~ C L)) where c^ a are the annihilation operators of the electrons in 
the lead and t is the tunnelling amplitude. We assume a constant density of states po in 
the lead such that the (bare) tunnelling probability is given by T = 2nt 2 po. We have 
calculated the differential Andreev conductance G(V) using exact diagonalization for 
an open chain of length N = 24 taking only the two lowest energy states (with different 

not yet been fully established |22) . 




Figure 3. Tunnelling conductance along the paths (a) c\ and (b) ci shown in figure [5] 
The arrows indicate the approximate phase boundaries. 



fermion parity) into account^ In figure ^ a), the broadened conductance 

J-Aw/2 Aw 

is plotted along c\\ the broadening Au ~ max(eV, ksT) is given by the maximum of 
bias voltage V and temperature T. For strong coupling T to the lead as well as small 
broadening Au, the phase boundaries are clearly visible as points where the conductance 
jumps from one to zero and vice versa. Weaker coupling to the lead will lead to 
a suppression of the conductance in the topological phase, i.e., unitary conductance 
cannot be observed. On the other hand, a larger broadening will eventually smear 
out all transitions. Thus in order to enable an experimental detection of the phase 
boundaries the energy broadening should be small (in units of Em) while the coupling 
to the lead has to be sufficiently strong. Recent experiments on proximity coupled 
nanowires [4j indicate that a Majorana coupling Em /kg — 100 mK is realistic, which is 
well in the range of experimental accessible temperatures. The Majorana coupling sets 
the topological gap A T (see inset of figure [2]) as the other couplings (Ej and E c ) can 
be designed in a large parameter range |21j . 

In order to study the regime of strong interactions we use the set-up with individual 
back gates and induced charges = (— l) l q. In this set-up we consider the path 
c 2(q) = Q-V, Ta cos{nq/e)) with IV = 0.3 Ta = 1.5 Em (see the green curve in figure EJ). 
The conductance G along c<i is shown in figure [3]^b). At q = the path starts deep in the 
band insulator phase and enters the topological phase at q ~ 0.2 where the conductance 
becomes nonzero. When further increasing q we observe weak oscillations which are 

% The two-level approximation is appropriate in the topological phase (where we have two levels 
separated by 3> S from the rest), as well as in the trivial phase (where there is a unique ground 
state) and in the Mott phase (where there are two degenerate ground state with the same fermion 
parity) as S > Auj in the latter cases such that the current vanishes. 
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due to the finite size lifting 5 of the ground state degeneracy. For larger values of q the 
conductance stays zero in the incommensurate and commensurate CDW phases with 
the exception of the phase transition between them, where the conductance is nonzero 
due to finite-size effects. 

Above we have shown how to realize the ID ANNNI model using Josephson junction 
arrays and how to map out its phase diagram by measuring the tunnelling conductance. 
In this sense our proposed set-up constitutes a quantum simulator for the ID ANNNI 
model. In particular, the experimental control over the system parameters like the gate 
voltages allows to study the effects of disorder on the phase diagram. 

6. Relation to nanowires 

Interacting electrons in proximity-coupled semiconducting nanowires are described by 
the microscopic Hamiltonian [8] 



with iff(x) = fiff(x), ^i(x)) T the electron field operator, m the electron mass, \x the 
chemical potential, a the strength of the spin-orbit coupling, Ez the Zeeman energy 
due to the applied magnetic field, A s the s-wave pairing amplitude, and U the (short- 
range) Coulomb interaction. For sufficiently strong Zeeman energy (compared to the 
other energy scales), we only need to consider a single band similar to the ANNNI model 
discussed above. However, projecting the Hamiltonian (J7J) onto a single band strongly 
reduces the effect of the interaction. Specifically, we find for the effective interaction 
strength in the single-band model U/Em = mU^a 2 jf?LE\ <C 1 with L the length of 
the nanowire. In this way, interacting nanowires subject to a strong Zeeman field are 
always in the weak coupling regime [5]. In contrast as we showed above, the strong 
interaction regime for spinless fermions is readily accessible in the case of nanowires in 
Josephson junction arrays. 

7. Conclusions 

We have analysed a ID array of Josephson junctions featuring Majorana modes, where 
capacitances between adjacent islands lead to interactions between the Majorana modes. 
We have shown that repulsive interactions generically facilitate the topological phase 
due to their competition with the on-site charging energies. Finally, we have proposed 
a tunnelling experiment to detect the phase boundaries. 
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Appendix A. Derivation of the quantum phase slip rate 

In this appendix, we present more information about the derivation of the rates 
and r^. Starting with Ta, we are interested in the event that the phase 4> k of a 
single island changes by 2tt. The relevant tunnelling matrix element t& oc (27r|e~* m |0) 
can be evaluated in the path-integral formalism by going to imaginary time r = it, 
cf. [231 El E5], 

tAOC [ V\ct> k }e- h '^ dTLE (A.l) 



subject to the boundary condition 0^.(0) = and A(f) k = ^fc(oo) — 0^(0) e 27rZ\{0}. To 
exponential accuracy, the path-integral is dominated by the classical paths <j) c k l ' n which 
minimize the action Se = / dr Le, i.e., 

t A ~j2 e ~ SE[4>t ' n] ( A - 2 ) 

n 

where n is an index enumerating the different paths in the case that there are different 
minima of the action. 

For i] = 0, the relevant part of the action is well- approximated by (keeping only 
terms which depend on 

T Vj 
oo ' ' v 



SeM = I dr - % ±4f k + ^(1 - cos0 fc ) , (A.3) 



o 



L E 



where ' denotes the derivative with respect to r and we have neglected the potential 
proportional to E M Ej. As the action does not depend directly on r, the energy 
along the classical path minimizing the action is conserved, 

E = -L E = T-Vj. (A.4) 

For r = we have T — Vj — such that E = in our case. 

We can express the kinetic energy in terms of the potential and obtain 

T= ik^ = E + Vj (A ' 5) 

with which we can get an alternative expression for the measure (the capacitance matrix 
acts as a metric) 



dT = ^ 16E c ( h E + V J ) ^ (A ' 6) 
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Due to the conservation of energy, we can go over to the Euler-Maupertuis action 
So = Se + JdrE (note that in our case E — such that So is i n fact equal to Se) which 
can be rewritten employing ( 1A.6I) as 



So[<Pk] = \l^j d ^ VETVj - ^A0 fc (A.7) 

which is independent on imaginary time and only depends on the path chosen [26J. 

The action is minimized for A<p k = ±27r as each additional phase slip by 2n increases 
Sq. The expression corresponding to the first term in flA.7j) is independent on ±. The 
classical path corresponds to increasing (fi k by 2n such that, cf [T7] . 



J d(j> k y/E + Vj 



;i-cos0 fc ) = J— -. (A.8) 



4-E'c J V 4£ c Jo V E c 

We obtain the final result 



tA ~ e -V^j/^C e -^A0 fe /2e _ e -y/8Ej/E c cos(g7r/e ). ( A .g) 

A</. fe =±27r 

valid up to exponential accuracy. In the main text, we use the result t& = Ta cos(7rg/e) 
with Ta — E x J,^E i J^e~^ %Ej l Ec . In fact the charge q should be replace by q + en^. The 
reason is that due to the Majorana term the action is not 2n but only 47r periodic in <frk- 
In the calculation above, we however assume the action to be 2ti periodic. In fact, the 
action can be made 2ir periodic by a gauge transformation on the expense of replacing 
q t— > q + en k . More information on this rather subtle point can be found in [T5], ITi] . 

The prefactor E^Ej of depends on the shape of the potential close to the 
turning points <fi k ~ 0, ±2n and cannot be obtain in our simple semiclassical analysis 
which only captures the physics up to exponential accuracy. However, the scaling of 
the prefactor can be obtained from summing up the instanton contributions [21] or by 
matching it to the exact solution of the Mathieu equation \T7\. In our case, the potential 
always is given by Vj ~ \Ej<\>\ for <\> k <C 1 such that the same prefactor EqEj (from 
the Mathieu equation) appears for all the tunnelling amplitudes. 

In the case r] ^ 0, it is important to notice that the phases on the different islands 
do not completely decouple. In lowest order in 77, we need to take the phases 4> k ±i on 
the islands k ± 1 into account. Due to the symmetry of the problem, we know that 
0fc-i( r ) = <Pk+i{ T )- The relevant part of the action reads 

Sk= I dr ^(4 2 + 2(l-rM 2 +1 ) 



16E C 

- ^(0' fc + 2<P' k+l ) + ^(3 - cos0 fc - 2cos0 fc+1 )] . (A.10) 

Vj 

Following the same line of calculation as going from (1A.3j) to (1A.7j) . we obtain 



So = ^AE^J d<Pk V [1 + 2(1 _ 'M+iC^) 2 ]^ - J ^ ) ( A - 1 1 » 



•Ccff 
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where we expressed the path by giving fc+1 as a function of <p k . As we are interested 
in processes where <p k changes by A<ft k = ±2ir, we need to find <fik+i {4>k)'- [0, ±27r] i— > M. 
with 0fc+i(O) = 0fc+i(±27r) such that the action is minimized. We will find the solution 
which corresponds to A(p k = 2ir below. The second solution with A<p k = — 2tt can by 
obtained via the symmetry <pi y — 0j, Vi of the Lagrangian. 

As the second term in flA.101) is independent of the path (it depends only on the 
boundary condition), we only need to minimize the first term. The extremum is attained 
when the function 4>k+i{4>k) fulfils the Euler-Lagrange equation 
d dC c s dC c s 



d(j) k d(p' k+l d(f) k ' 

To first order in r] [assuming G O 

4(1 - cos0 fc )0' fc ' +1 + 2 sin 

Employing the substitution (p k+ i((f) k ) 
order equation in /' of the form 

2cos(0 fc /4)[8sin 2 (0 fc /4)r(0 fc ) - r/cos(0 fc /2)] 

+2[3sin(0 fc /4) + sin(30 fc /4)]f (0 fc ) = 

which can be integrated with the solution 

logcos(|(/> fc ) 



(A.12) 

, the equation assumes the form 

*) 0Ui ~ V sin ~ 20 fc +i = 0. 
tan(0fc/4)/(0fc) reduces the equation to first 



f{<j) k ) = 77 log 



- sm(-(b k ) 
2 v 2 ' 



-71- 



sin 2 (i0 fc ) 



(A.13) 



(A.14) 



Plugging the solution into (lA.lljl and retaining the first nonvanishing term in r\ yields 



So = \/8Ej/E, 



c 



1 



7T" 



12 



96 



V 2 + 0( V 4 



iqA(f) k 



(A.15) 



2e 

Summing up the two contributions with A<p k = ±2tt, we obtain a term 
proportional to cos(7rg/e) as before, with the proportionality constant given by = 

E lfi £ 3/ V V 8 ^/ £ o[l+(T 2 -12)r? 2 /96] _ 
C J 

For 1] ytz 0, we get additionally a next-nearest neighbour interaction due to phase 
slips where both (fi k and (fi k+ i change by 2n. In fact, due to the symmetry of the problem, 
we can set (f) k +i(r) = <p k (T). The term of the action which change with <p k and <p k +i 
are given by ( 1A.3I) for each of the islands and the additional contribution of the cross 
capacitance. Thus, we have 

-m-v) ,, 2 



s f 



dr 



WE, 



c 



+ 2Ej(l-cos(j) k ) 



(A.16) 



which leads to 



So 



(2 - V )Ej 



2E, 



c 



1 — COS (j) k 



iq(A<f> k + Acf), 



k+lj 



2e 



4y/(2-rj)Ej/E, 



c 



iq(A(j) k + Acf) k+1 ] 
2e 



(A.17) 



As A(f) k+ i = Acj) k (because the two phases slip together), the tunnelling amplitude thus 
assumes the form IV cos(2vrg/e) with IV = E^Ej^e' ^™V-v)Ej/e _ 
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